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Die Bedeutung des Konkurrenzfaktors bei der 
Stabilisierung historischer Arealgrenzen 


G. PILS 


Abstract: Alpine botanical chorology in the last decads has been dominated by hi- 
storical aspects. Therefore recent area restrictions without evident ecological reasons — 
unfortunately a very frequent case — mainly are ascribed to „genetic impoverishment“ 
due to glacial area contractions (for more detailed critics see PILS 1988). As these histo- 
rical theories suffer from a chronical lack of falsification possibilities, this paper attempts 
to combine the recent advances in population biology with well established chorological 
facts. 


Relying e. g. on ENDLER (1977), the theoretical basis of our models is the assumption, 
that nowadays distribution boundaries are maintained by the combined effect of two op- 
posite forces: Gene flow by means of diaspore transport tends to expand distribution 
areas, whereas competition with already established species counteracts this tendency. 
These assumptions allow the formulation of differential equations. By solving them re- 
peatedly, the evolution of a set of demes („Zellen“) can be simulated for thousands of 
years. 


The basic model is based on the assumption of a positively frequency-correlated fitness 
function (PFCFF) for each of two competing species. In nature this should be realised in 
allogamous species the diaspore production of which is (sometimes drastically!) reduced 
when they are grown at very low population densities. Experimental proof for this as- 
sumption comes from extensive self-pollination experiments in the genus Festuca 
(Poaceae) by AUQUIER (1977). In its mathematical and also biological consequences our 
model corresponds to the „tension-zone“ models of modem population genetics, which 
explain stable steps between subspecies by fitness reductions due to heterozygote 
disadvantage. The most interesting aspect of all these models based on PFCFFs is that 
the final position of the step between competing (sub-)species depends on the point of 
first contact, and not on environmental factors. Therefore, historical area constellations 
may be „frozen“ and recent distribution areas sometimes outline Pleistocene refugia. 


Wherever PFCFFs exist, isolated populations need a certain minimum area to be stable 
against invasion by gene flow (seed dispersal) from competing species. Extinction due to 
reduction of the minimum area below a certain threshold may be extremely slow at the 
beginning, but accelerates in an autocatalytic way in its final phases. Establishment of 
new mutants with some degree of reproductive isolation, e. g. polyploids, within existing 
populations with a PFCFF virtually should be impossible if one would not like to have 
recourse to additional hypotheses as e. g. a very effective vegetative propagation or an 
otherwise very superior fitness of the mutant. For explaining the high percentage of neo- 
polyploids in the Alpine flora, the importance of founder effects during expansion peri- 
ods is stressed. As shown by simulations, far-range dispersal of a mutant diaspore into 
still competition free areas (e. g. after ice retreat) rapidly leads to the forming of big po- 
pulations by logistic growth. Given enough time, such founder populations afterwards 
should be numerous enough to be stable against invasion of the original diploid and 
themselves be able to conquer new areas after further ice retreat. 


Key words: Chorology, population biology, evolution. — Computer simulation, 
tension zone, competition, relic species. 
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Einleitung 


Die botanische Arealkunde im Alpenraum ist spätestens seit den eingehenden 
Untersuchungen MERXMÜLLERS (1952-54) ausgesprochen historisch orientiert. 
Ursache dafür ist die auffällige Bindung einer Reihe lokal verbreiteter Pflanzen- (und 
Kleintier-) Arten an eiszeitlich gering oder gar nicht vergletscherte Räume, ohne daß 
diese Arealbegrenzungen durch rezente Ökologische Gradienten hinreichend zu 
erklären wären. Mit MERXMULLER! wird daher heute vielfach davon ausgegangen, 
daß es sich um „Relikte“ handelt, d. h. einst wesentlich weiter verbreitete Sippen, 
welche im Gefolge der eiszeitlichen Gletschervorstöße in weiten Teilen ihres 
einstigen Areals ausgelöscht wurden und nur in eiszeitlich wenig vergletscherten 
Refugialgebieten („massifs de refuge“) den Sprung in die Gegenwart geschafft 
hätten. Die bemerkenswerte Fixierung der vom Gletscher aufgezwungenen 
Arealgrenzen bis heute, d. h. etwa 15.000 Jahre nach dem Gletscherrückzug, wird 
dabei meist dadurch erklärt, daß sich manche dieser Sippen, an Biotypen verarmt, als 
nicht mehr genügend wanderungsfähig erwiesen, sodaß ihre glazial gestalteten Areale 
bis heute konserviert blieben. 


Ein Hauptkritikpunkt gegen dieses historizistische Modell der Arealformung ist seine 
mangelnde Falsifizierbarkeit (PILS 1988, 1994). Einerseits besteht keine realistische 
Aussicht, die oft nur postulierte, im besten Fall aus aktuellen Arealbildern 
erschlossene, einstige weite Verbreitung der als „Relikte“ apostrophierten Sippen 
tatsächlich durch Fossilbelege nachzuweisen. Andererseits wurde bisher allenfalls in 
Ausnahmefällen die oft strapazierte Biotypenverarmung als Grundlage der aktuellen 
Ausbreitungsschwäche von Reliktarten tatsächlich durch konkrete Untersuchungen 
belegt (z. B. bei den extrem individuenarmen, voneinander weit isolierten und sich 
noch dazu hauptsächlich vegetativ vermehrenden alpischen Reliktpopulationen von 
Saxifraga cernua L., vgl. BAUERT 1994). 


Als ein möglicher Ausweg aus dieser Sackgasse bietet sich die computergestützte 
Simulation populationsbiologischer Modelle an, wie sie etwa bereits vielfach und mit 
überzeugenden Ergebnissen zur Klärung analoger Probleme in der Populationsgene- 
tik eingesetzt werden (z. B. ENDLER 1977). Ziel der vorliegenden Arbeit soll es daher 
sein, derartige Modelle auch zur Überbrückung des virulent gewordenen Gegensatzes 
zwischen historischer und aktualistischer Sicht der Arealbildung bei Gefäßpflanzen 
einzusetzen. 


I MERXMÜLLER (1954: 132, Punkt 8 der abschließenden Zusammenfassung): „Rezente Arealbe- 
schränkungen auf einzelne Erhaltungszentren (wie das Nordareal der Arten mit Süd-Nordost- 
Disjunktion) wie auch Disjunktionen innerhalb der Nordalpen (so die „Bayerische Lücke“) werden 
durch mangelhafte Rückwanderung und reduzierte Ausbreitungstendenz erklärt. Nur Arten mit 
überdurchschnittlichem Ausbreitungsvermögen nehmen heute wieder zusammenhängende Areale 
ein“. 
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Material und Methode 


Die Simulationen wurden mit einem DOS-kompatiblen PC unter Zuhilfenahme des 
handelsüblichen Tabellenkalkulationsprogramms Excel 5.0 durchgeführt. Ein Ar- 
beitsspeicher von mindestens 4 Mb erwies sich dabei von Vorteil. 


Theoretische Grundlagen 


Basis der vorliegenden Untersuchungen ist die Hypothese, daß es sich bei sämtlichen 
derzeit schlecht verstandenen und daher ins nicht falsifizierbare Feld der Vergangen- 
heit abgeschobenen ,,reliktischen“ Arealgrenzen letztlich um konkurrenzbedingte 
Phänomene handelt. Folgende Hinweise liegen dazu bereits vor: 


l. Wie bereits in PıLs (1988) an Hand konkreter Beispiele erläutert, fehlt eine ganze 
Reihe weitverbreiteter Alpenpflanzen gerade den Refugialräumen. Dement- 
sprechend liegen bisher keinerlei konkrete Nachweise dafür vor, daß sich Refugial- 
räume und vegletscherte Gebiete in ihren absoluten Artenzahlen pro Einheitsfläche 
signifikant voneinander unterscheiden2. Dies deutet auf einen allgemeinen und 
letztlich konkurrenzbedingten Vikariismus zwischen Refugial- und Nicht-Refu- 
gialarten hin. Der historizistische Ansatz eignet sich zur Erklärung dieses Phäno- 
mens nur schlecht: Nimmt man nämlich vor den eiszeitlichen Verwüstungswellen 
in allen Alpenteilen etwa gleiche Artenzahlen pro Einheitsfläche an, so wäre heute 
in den „eiszeitlich gereinigten“ zentralen Alpenteilen eine spürbar artenärmere 
Flora zu erwarten als in den Refugialgebieten?. 


2. Auch die hauptsächlich von zoologischer Seite bereits durchgeführten mathema- 
tisch-experimentellen Analysen von Übergangszonen zwischen hybridisierenden 
(Unter-) Arten beruhen sämtliche auf dem Grundpostulat, daß ein Ausbleiben der 
in diesen Fällen auf Grund des Genflusses zu erwartenden Vermischung auch über 
lange Zeiträume nur durch gegengerichtete Selektionsvorgänge verhindert werden 
kann (vgl. zuletzt etwa HARRISON 1993). Das einzige Hindernis für die Arealaus- 
weitung einer Subspecies ist also auch hier die Präsenz einer konkurrierenden an- 
deren Unterart. 


2 Derartige Untersuchungen wurden zumindest im Alpenraum noch nie systematisch durchgeführt. 
Zu beachten wäre dabei jedenfalls die oft sehr unterschiedliche Vertikalerstreckung der zu verglei- 
chenden Gebirgsmassive. So dürfte etwa der oft gerühmte Artenreichtum des Alpensüdrandes wei- 
testgehend auf eine Höhenstufenerweiterung nach unten zurückgehen, welche das Auftreten zahlrei- 
cher submediterraner Elemente ermöglicht. 


3 In der Praxis hilft man sich hier regelmäßig mit Nunatakkem, also über die Eismassen aufragen- 
den Bergen, welche auch in stark vergletscherten Alpenteilen ein Überdauern artenreicher Gesell- 
schaften bzw. „reliktischer‘‘ Arten ermöglicht hätten. In der Praxis ist aber selbst beim Vergleich 
mit den modernsten heute zur Verfügung stehenden Karten der eiszeitlichen Vergletscherung (etwa 
Husen 1987) eine Falsifizierung solcher hypothetischer ones (prinzipiell!) niemals 
möglich. Sogar bei einer völlig fehlenden Korrelation zwischen aktueller Pflanzenverbreitung und 
würmzeitlich unvergletschert gebliebenen potentiellen Nunatakkerstandorten läßt sich nämlich den- 
noch mit dem Postulat einer „Rückwanderung auf kurze (mittlere, weite) Distanz“ jede beliebige 
Arealkonfiguration historisch erklären. 
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Nach der vermuteten Ursache dieser Selektionsprozesse lassen sich die bisher prä- 
sentierten Modelle in zwei Hauptgruppen unterteilen: 


1, Selektions-Gradienten-Modell („geographical-selection-gradient model“): Entlang 
eines räumlichen (,,exogenen“) Selektionsgradienten ändert sich die relative Fit- 
ness der Genotypen (Unterarten). Die Folge ist die Ausbildung einer räumlich be- 
grenzten und zeitlich stabilen Übergangszone (Hybridzone, „step“) zwischen den 
Genotypen. Dieser Ansatz wurde bereits von FISHER (1950) mathematisch formu- 
liert und ın der Folge insbesondere von SLATKIN (1973, 1975), MAY & al. (1975), 
sowie ENDLER (1977) mittels Computersimulationen eingehend ausgetestet. Ob- 
wohl dieses Modell das von vornherein einsichtigste ist, krankt doch seine prak- 
tische Umsetzung sehr häufig daran, daß insbesondere in den lokalklimatisch und 
auch edaphisch extrem kleinräumig strukturierten Hochgebirgen die Suche nach 
überregionalen Umweltgradienten, welche für großräumige Fitnessunterschiede 
verantwortlich gemacht werden könnten, oft ergebnislos bleibt. 


2. „Lension zone model“ sensu BARTON & HEWITT (1985), basierend auf negativ fre- 
quenzkorrelierter Selektion 


a) Hybridunverträglichkeitsmodell (z. B. BARTON 1979a,b, BARTON & HEWITT 
1989): Primäre Ursache für die Ausbildung enger und konstanter Übergangszonen 
ist in diesem Fall die („endogene“) Selektion gegen intermediäre (= Hybrid-) Ge- 
notypen auf Grund deren reduzierter Fitness (,,heterozygote disadvantage‘). 


b) Genfrequenz-Selektionsmodell (,,frequency-dependent selection“, MALLET 1986, 
1989, 1993): In diesem Fall richtet sich die Selektion aus „exogenen‘“ Gründen ge- 
gen das jeweils seltenere von zwei konkurrierenden Allelen bzw. Genkombina- 
tionen. Bisher wurde dieser Fall nur bei südamerikanischen Heliconius-Arten 
(Lepidoptera) beschrieben. Hier liegen Rassenkreise von durch Müllersche Mimi- 
kry geschützten ungenießbaren Arten vor, bei denen seltene Varianten potentiellen 
Freßfeinden noch nicht bekannt sind. Sie werden daher - trotz ihrerer Ungenieß- 
barkeit — häufiger attackiert als verbreitete Formen. Auf diese Art wird natürlich 
ebenfalls die Ausbreitung einer Rasse ins Areal einer anderen verhindert. 


Obwohl sich in ihrer mathematischen Behandlung alle 3 Ansätze sehr ähnlich sind, 
bestehen in ihren arealkundlichen Konsequenzen grundlegende Unterschiede zwi- 
schen 1 und 2. Während bei den Selektions-Gradienten-Modellen letztlich die Um- 
welt die Lage der Arealgrenze zwischen einem vikariierenden Sippenpaar diktiert, 
und diese sich nach Störungen daher stets wieder auf ihrem alten Wert einpendelt 
(permanente Umweltveränderungen natürlich ausgeschlossen), ist dies bei „tension 
zones“ nicht der Fall. Wird die aktuelle Arealgrenze verschoben, was etwa durch eis- 
zeitliche Arealverwüstungen in großem Maßstab passiert sein muß, so pendelt sich 
nach dem Ausklingen der ursprünglichen Störung die neue Lage der Hybridzone am 
ersten Kontakt nach erfolgter Rückwanderung ein. 


Insgesamt befassen sich alle bisher erwähnten Modelle mit reproduktiv noch nicht 
isolierten und daher meist als Subspecies eingestuften Sıppen. Kaum Ziel konkreter 
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Untersuchungen waren dagegen bisher Arealgrenzphänomene zwischen (oft heftig!) 
konkurrierenden, aber reproduktiv isolierten Arten. Dabei liegen auch hier durchaus 
konkrete Hinweise für die entscheidende Bedeutung der Konkurrenz bei der Stabili- 
sierung von Arealgrenzen vor. 


Als ein exemplarisches Beispiel aus den eigenen Untersuchungen seien etwa die 
beiden Violettschwingelsippen Festuca picturata PILS und Festuca puccinellii PARL. 
erwähnt, deren Verbreitungsgebiete sich — abgesehen von einer auffällig schmalen 
Überlappungszone im Bereich der Zillertaler Alpen — komplementär verhalten (PILS 
1980). Beide Arten haben nach dem bisherigen Kentnisstand weitgehend identische 
Ökotopansprüche und kommen innerhalb ihrer Areale häufig, ja oft sogar 
bestandbildend vor (Trifolio-Festucetum violaceae BR.-BL. 26 und verwandte 
Gesellschaften). Typische Standorte sind tiefergründige Hangmulden in der (sub-) 
alpinen Stufe über neutralen bis sauren Böden (,,Urwiesen“). Hybriden zwischen den 
beiden auf verschiedener Ploidiestufe stehenden Sippen (Pseudovikarianten sensu 
FAVARGER 1962) sind bisher nicht bekannt geworden (vgl. auch KIEM 1987). Auch 
die bisherigen morphologischen Befunde sprechen gegen eine besonders enge 
Verwandtschaft® . ’ 


Bereits in früheren Arbeiten äußerten wir die Vermutung, daß die Ursache für die 
derzeitige Arealgrenzlinie in einer konkurrenzbedingten Blockade als Resultat einer 
gegensinnigen Rückwanderung nach dem Ende der letzten Eiszeit zu suchen sei (PILS 
1980, 1988). Ähnliches gilt für eine Reihe vikariierender (Unter-)Arten aus dem Tier- 
reich bereits als hinreichend belegt, wobei intensive populationsbiologische Unter- 
suchungen bereits eine Fülle hochinteressanter Zusammenhänge ans Licht gebracht 
haben. Zu den bekanntesten dieser auch Österreich tangierenden Sippenpaare gehö- 
ren Raben- und Nebelkrähe (Corvus corone corone L. - C. c. cornix L., z. B. 
AUBRECHT 1979) sowie Rot- und Gelbbauchunke (Bombina bombina L. — B. variega- 
ta L., z. B. SZYMURA 1993). Nachfolgend soll versucht werden, ein „tension-zone 
model“ auch für einen derartigen, letztlich zu geographischer Vikarianz führenden 
Konkurrenzausschluß zwischen allogamen, nicht hybridisierenden Pflanzenarten zu 
erarbeiten. 


Allogamie als Grundlage positiv frequenzkorrelierter Fitness bei Blütenpflanzen 


Viele Pflanzenarten, und zwar insbesondere solche längerfristig stabiler 
(nichtruderaler) Standorte, haben markante Selbstbestäubungshemmnisse entwickelt. 
Im Rahmen unserer eigenen biosystematischen Untersuchungen an der Gräsergattung 
Festuca L. konnten wir reduzierten Samenansatz an ısolierten Kulturexemplaren wie- 


4 Die Ursprünge der hexaploiden Festuca puccinellii sind aller Wahrscheinlichkeit nach bei der 
bisher ausschließlich als diploid bekannten Festuca violacea SCHLEICHER ex GAUDIN der Westalpen 
zu suchen, während die diploide E-alpisch-karpatische Festuca picturata große Affinitäten zur 
ebenfalls diploiden, aber streng calcicolen SE-alpisch-südkarpatischen Festuca nitida KIT. ex 
SCHULTES aufweist. 
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derholt feststellen. Ausführlich statistisch belegt wurde dieses Phänomen bei Festuca 
aber schon von AUQUIER (1977, Tab. 1). 


Wie aus beigefügter Tabelle hervorgeht, sind zumindest bei fast allen daraufhin bis- 
her eingehender untersuchten Festuca-Arten isolierte und damit zur Autogamie ge- 
zwungene Einzelpflanzen ın ihrer Fruchtbarkeit und damit ihrer relativen Fitness im 
Vergleich zu Pflanzen im Bestand extrem eingeschränkt. Dies gehi soweit, daß schon 
AUQUIER (1977) anmerkte: „Anemophile Populationen brauchen tatsächlich eine 
gewisse Pollenmenge, um überleben zu können. Mit einigen Tests haben wir im Ex- 
perimentalgarten beobachten können, daß eine Entfernung von einigen Metern schon 
die Kreuzbefruchtung bei den Festucen verhinderte“. 

Nur bei den beiden untersuchten Waldarten Festuca gigantea und Festuca altissima 
tritt dieses Phänomen praktisch nicht auf; wie bereits von AUQUIER vermutet wohl 
deshalb, weil beide Arten auf Grund der durchschnittlich geringeren Windstärken an 
ihrem Naturstandort schon von Natur aus stärker auf Autogamie angewiesen sein 
dürften. 
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| 0,68 (1,22) 26,87 (13,62) 0,96 (0,07) 
4,16 (5,54) 34,99 (13,47) 0,77 (0,28) 
0,62 (0,70) 42,69 (17,94) 0,97 (0,03) 


Festuca stricta 6,01 (6,64) 43,94 (18,75) 0,78 (0,21) 
HOST. ssp. trachy- 

phylla (HACK.) 

PATZKE 

Festuca nigrescens 3,11 (4,44) 42,63 (21,14) 0,88 (0,15) 
LAM. 

Festuca pratensis 3,55 (2,51) 53,01 (18,44) 0,88 (0,05) 
HUDS. 

Festuca arundinacea 2,85 (1,62) 36,82 (10,46) 0,86 (0,08) 
SCHREB. SSp. arun- 

dinacea 

Festuca heterophylla 0,90 (1,56) 28,33 (5,39) 0,95 (0,09) 
LAM. 

Festuca gigantea 83,00 (0,71) 89,20 (3,68) 0,04 (0,02) 
(L.) VILL. 

Festuca altissima 86,8 (eine Population) ee fF O 
ALL. 


Tab i: Berechnung relativer Fitnessunterschiede (d) aus den von AUQUIER (1977) ermittelten 
Karyopsenerträgen isolierter bzw. fremdbestäubter Individuen verschiedener Festuca-Arten. Nach 
(3) liegen die Werte von d im Bereich 0<d<1. Ein d-Wert von 1 ergibt sich, wenn eine zur Auto- 
gamie gezwungene Einzelpflanze keinerlei Karyopsen erzeugt (entsprechend einer Fitness von 0). 
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Umgesetzt in die Modelle der Populationsdynamik, ergibt sich aus obigem Zusam- 
menhang zwingend, daß von zwei nischengleichen und dadurch maximal konkurrie- 
rende Sippen jeweils die seltenere auf Grund der geringeren Fremdbestäubungswahr- 
scheinlichkeit in ihrer Fitness im Vergleich zum häufigeren Konkurrenten zusehends 
ins Hintertreffen gerät, und zwar umso mehr, je seltener sie bereits geworden ist. Die 
Situation entspricht damit prinzipiell der einer populationsgenetischen „tension zone“ 
sensu BARTON & HEWITT (1985). Damit muß die für unsere Überlegungen wichtigste 
Folge eines derartigen autokatalytischen Prozesses jedenfalls sein, daß sich zwei in 
sekundären Kontakt geratene Sippen (etwa nach postglazialen Expansionen) nicht 
vermischen, sondern im Gegenteil unter der Bildung von vergleichsweise scharfen 
Übergangszonen arealmäßig weitgehend getrennt koexistieren sollten. 


Allgemeine Grundhypothese 


Gelangen zumindest teilweise selbststerile, nischengleiche (= total konkurrierende) 
Pflanzenarten z. B. durch postglaziale Rückwanderung in sekundären Kontakt, be- 
ginnen zwei gegenläufige Kräfte zu wirken: Einerseits tendiert kontinuierlicher Sa- 
mentransport (= Genfluß) dazu, beide Nachbarn ineinander wandern zu lassen. Ande- 
rerseits eliminiert die negativ frequenzkorrelierte Selektion ständig selektiv Individu- 
en der jeweils selteneren Art und wirkt damit der Vermischung entgegen. Letztlich 
wird sich an der Arealgrenze ein pseudostationärer Gleichgewichtszustand einstellen, 
d. h. eine mehr oder weniger enge Mischzone, ın der die entmischende Kraft der Se- 
lektion und die diffusionsartig mischende Gegenkraft des Genflusses einander aufhe- 
ben. Diese Bewegung zum Gleichgewicht hin läßt sich folgendermaßen mathema- 
tisch formulieren: 


Modellannahmen für die Computersimulation 


e Zwei reproduktiv völlig isolierte Arten A und B befinden sich in (sekundärem) 
scharfem Kontakt. 


e Innerhalb des Gesamtareals beider Arten treten keine umweltbedingten Fitness- 
unterschiede zwischen beiden Gegenspielern auf (entsprechend den bei nicht re- 
produktiv isolierten Sippen angewandten „tension-zone models“). 


e Die Gesamtpopulation setzt sich aus linear gleichmäßig angeordneten, gleichstar- 
ken Teilpopulationen (Zellen, ,,demes“) zusammen, zwischen denen ein wechsel- 
seitiger Genfluß durch Samentransport besteht (= „linear stepping-stone model“). 


e [A], [B] = Frequenz (Prozentanteile) von A, B in der gleich 100 gesetzten Zelle 
e Innerhalb dieser Zellen stehen A und B in totaler Konkurrenz: 
[3] = 100 - [A] (1) 
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Fitnessfunktionen (Abb. 1): 
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Abb. 1: Fitness von A und B (W, und W,) bei Annahme eines linearen, positiv frequenzkorrelierten 
Zusammenhanges. 


Fitness von A: W = (1-d) + (d/50*[A]) (2a) 
Fitness von B: W,= (1-d) + (d/50*[B]) (2b) 
oder nach (1) W,= (1-d) + (d/50*(100-[4])) . (2c) 


d = maximaler relativer Fitnessunterschied zwischen beiden Sippen: 0<d<] 


In der Praxis berechnet sich d folgendermaßen aus experimentell gewonnen Fitness- 
messungen (z. B. dem Samenansatz selbst- bzw. fremdbestäubter Individuen): 


d=1-( Wan! (Wain Wna)!2)) (3) 
W „= Minimalfitness, W_ = Maximalfitness 


Für praktische Untersuchungen wären aber jedenfalls die durch Selbstbestäubung 
entstandenen Diasporen zusätzlich noch in Hinblick auf eine etwaige verminderte 
Keimfähigkeit oder gesenkte Vitalität der daraus entstandenen Individuen zu unter- 
suchen. Bemerkenswert sind in diesem Zusammenhang auch die Beobachtungen 
CONNORS (1957) an Festuca novae-zelandiae (HACKEL) COCKAYNE, wonach durch 
Selbstbestäubung entstandene Individuen untereinander gekreuzt keine Nachkom- 
menschaft mehr ergeben, mit anderen Kreuzungspartnern aber sehr wohl. Erste dies- 
bezügliche Versuche von AUQUIER (1977) mit Festuca filiformis POURRET, F. le- 
manii BAST. (= F. bastardii KERGUELEN & PLONKA) und F. stricta HOST ssp. trachy- 
phylla (HACK.) PATZKE (= F. brevipila TRACEY) bestätigten die Angaben CONNORS. 


Ableitung der Simulationsformeln 


Ohne Berücksichtigung des Samenaustausches mit den Nachbarzellen errechnen sich 
die relativen Frequenzanteile von A und B im Jahr n+1 folgendermaßen: 


[4,..] #,,°[4,] (4a) 


[B,..] = W,,*[B,] (4b) 
W, = Fitness von A im Jahr n; W, = Fitness von B im Jahr n 


Im Grundmodell beschränkt sich der Genfluß mit der Umgebung auf rechte und linke 
Nachbarzelle. 
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[4,1] = (1-28)*[4,]* Wun Fe" AL We PS TA Wan (5a) 
(B,..) = (-28Y {8,1 Wn 8B], WM nn (Sb) 
g = Genflußfaktor (Anteil der Neupopulation, welcher durch Samenanflug aus der 
betreffenden Nachbarpopulation stammt) 
[4,],= Frequenz von A in der rechten Nachbarzelle zum Zeitpunkt n 
[4,],= Frequenz von A in der linken Nachbarzelle zum Zeitpunkt n 


Unter Berücksichtigung von (1-4) stehen daher [4] und [B] zum Zeitpunkt n+1 in fol- 
gendem neuen Verhältnis: 


[4,1 : (B,.,] = (-2g)*((1-d) + (d/50"[An]))*14,]}+ g*[4,],*((1-d) + (d/50*[A,],))+2* 
[4,1], *((1-d)+(4/50°(4,1)):(1-2g)*((1+4)+(d/50*(100-[4,7)))*(100-[4,1+8*(100-[4,],)* 
((itd)+ (d/50*(100-[4,],))) +g*(100-[4,1,*((1+4)+(4/50*(100-[4,],))) 

Wobei nach (1) erneut gelten muß: [4,,,]+[B8.,,] =100 


Daraus läßt sich sich folgende allgemeine Formel für die Berechnung von [A 
ten! 


[A4, = 100/((((((100-[4,])*(1-2g)*((1-d)+(d/50*(100-[4,}))))+(g*(100-[4,])*((1- 
d)+(d/50*(100-[4,],))))+(g*(100-[4,],)*((1-d)+(d/50*(100-[4,],)))) /(((1-2g)*[4,] *((1- 
d)+(d/50*[A,})))* @*[A,],* (A-d) +(a/50*[4,],))) +(e" 14) "(C1 -d) m 


] ablei- 


atl 


Fitnessunterschied (d) und Weite der Uberlappungszone 
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Abb. 2: Weite der Uberlappungszone bei verschiedenen Fitnessunterschieden und konstantem Gen- 
fluß (g = 0,1).- Für d-Werte von 0,5 bis 0,005 wurde die Simulation bis zur Erreichung des Gleich- 
gewichtes durchgespielt, für d = 0,0005 und d = 0,0 bis zur Generation 10.000. Die Kurve für d = 
0,0005 ist auf Grund des zu kleinen Areals von B nicht stabil. Im verwendeten Maßstab sichtbar 
wird dies aber erst bei Generationszahlen >> 10000 (vgl. Abb 5). Für d = 0,0000 würde der Gleich- 
gewichtszustand einer waagrechten Linie entsprechen. Die Annäherung daran erfolgt asymptotisch. 
Der Kurve d = 0,005exp liegt die exponentiell sinkende Fitnesskurve von Abb. 3 zugrunde. 
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Die Mathematik sigmoidaler Clines, welche durch ein Gleichgewicht von Selektion 
und Dispersion aufrechterhalten werden, wurde von populationsgenetischer Seite be- 
reits eingehend ausgelotet. Wir beschränken uns hier daher auf die Illustration der 
wichtigsten praktischen Konsequenzen unseres Modells. 


Wird die Weite der Ubergangszone (w) definiert als 1/max. Kurvensteigung, so gilt in 
Anlehnung an ENDLER (1977) und BARTON & HEWITT (1989) folgende Proportiona- 
lität: 


2. 
w = }/k „= const. * EE (7) 


k,. = max. Kurvensteigung; x = durchschnittliche Verbreitungsdistanz pro Genera- 
tion, g = Wahrscheinlichkeit eines Samenimportes aus sämtlichen Nachbarzellen. 


In dieser ersten Simulationsserie (Abb 2) wurde zunächst die Größenordnung von d 
(= maximaler relativer Fitnessunterschied) schrittweise um je eine Zehnerpotenz ab- 
gesenkt. Die Weite der Übergangszone verhält sich in diesem Fall proportional zu 


1/Va. 
w = const. * 1/Va 


Daß eine stärkere Selektion (größeres d) zu engeren Übergangszonen führt, ist in den 
allgemeinen Modellgrundlagen bereits vorgegeben. Umgekehrt verlängert sich mit 
abnehmender Selektionsintensität die Zeit bis zur Einstellung des Gleichgewichts 
sehr stark, was darauf zurückgeht, daß auf Grund der jetzt breiteren Übergangszone 
wesentlich stärkere Populationsverschiebungen relativ zur hier gewählten Ausgangs- 
position (Breite der Übergangszone = 0) zu bewältigen sind. 


Als bemerkenswerter Extremfall wurde auch die Kurve für d = 0,0 in die Graphik 
aufgenommen. Sie entspricht dem Fall, daß zwei gleich konkurrenzstarke Arten mit 
identischen Nischenansprüchen in Kontakt kommen, ohne daß irgend eine Art dichte- 
abhängiger Selektion auftritt. Die Annäherung an das endgültige Gleichgewicht 
(w= ) erfolgt hier asymptotisch (vgl. Abb. 6). Bei langlebigen Arten ohne Fernver- 
breitungseinrichtungen, also etwa auch bei unseren alpinen Schwingeln, könnte diese 
diffusionsartige Vermischung derart langsam fortschreiten, daß Änderungen inner- 
halb der von uns bisher einigermaßen dokumentierten Zeiträume gar nicht mehr 
nachweisbar sind. Die Versuchung wäre dann groß, auch hier von „eingefrorenen‘ 
Arealgrenzen zu sprechen. 


Letztlich wurde unter d = 0,005exp auch eine Fitnesskurve getestet, bei der mit sin- 
kender Populationsfrequenz die Fitness exponentiell abnimmt (Abb. 3). 
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Abb. 3: Mit abnehmender Populationsfrequenz exponentiell sinkende Fitnessfunktion (8). 


Die hier angewandte Fitnesskurve gehorcht folgender Gleichung: 


W, = 0,995+((1-(1-([4]/100))*)/100) (8) 
für 0 <[4]< 100 


Diese Fitnesskurve dürfte den tatsächlichen Verhältnissen insofern näherkommen, als 
bei Populationsanteilen >50 % das Ausmaß der Fitnessminderung durch Ausbleiben 
der Fremdbestäubung noch sehr gering dürfte. Konkrete diesbezügliche Untersuchen 
fehlen allerdings noch völlig. Der entscheidende Unterschied zur linearen Fitness- 
kurven ist jedenfalls die deutlich niedrigere Selektionsintensität bei mittleren Popula- 
tionsanteilen, die Auswirkungen auf die Weite der Übergangszone sind aber im Ver- 
gleich zu Änderungen bei d oder g gering. 


Genflußfaktor (g) und Weite der Überlappungszone 


Wie aus der allgemeinen Formel für die Weite der Übergangszone (7) abzuleiten, 
wirkt sich eine Erhöhung von x wesentlich stärker auf die Kurvensteigung und damit 
die Breite der Übergangszone aus als eine Steigerung von g. 


Bei konstanter Selektionsstärke (d = 0,005) wurden in Abb. 4 folgende Kombi- 
nationen von g und x simuliert: 


8,7 0,1, Erias 0,025 Ara 





Dabei wurde bei der Simulation einer höheren durchschnittlichen Verbreitungsdi- 
stanz (x) folgendermaßen verfahren: 


Kurve g,,: Drei Viertel der verbreiteten Samen wandert in die rechts und links direkt 
angrenzenden Zellen (g, = 0,15), ein Viertel aber jeweils in die beiden übernächsten 


(g, = 0,05). 
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Kurve g,,,,,: Die Hälfte der verbreiteten Samen wandem in die beiden angrenzenden 
Zellen (g, = 0,1), die andere Hälfte wandert zu gleichen Teilen in die nächsten 
4 Nachbarzellen (jeweils rechts und links; g,, g,, g,, g,jeweils 0,025). Die starke Ab- 
hängigkeit der Breite der Übergangszone von der Weite des Diasporentransportes 
wird bei dieser Annahme besonders augenscheinlich. 
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Abb. 4: Auswirkung der Variation des Genflusses auf die Weite der Uberlappungszone 
(fur d = 0,005). Nähere Erklärung im Text. 


Wird allgemein der Genfluß im Verhältnis zur Selektionsintensität zu stark gesteigert 
(oder umgekehrt natürlich die Selektionsintensität zu stark gesenkt, vgl. Kurve d = 
0,0005 in Abb, 4), so tritt der Fall ein, daß auf Grund der räumlichen Begrenztheit 
des Simulationsraumes die jeweils von vornherein weniger weit verbreitete Art von 
der arealmäßig überlegenen überrollt wird und dadurch aus dem Simulationsraum 
völlig verschwindet (Abb. 5). Dieser Vorgang hat autokatalytischen Charakter, d. h. 
er beginnt sehr langsam und schreitet mit zunehmender Unterwanderung der unterlie- 
genden Art immer rascher fort. Da solche Verdrängungsvorgänge an Barrieren mit 
verringertem Genfluß zum Stillstand kommen, wäre in der Praxis also zu erwarten, 
daß im limitierten Raum eines Gebirgsstockes sich jeweils nur eine von zwei über 
frequenzkorrelierte Fitnesskurven verbundenen Arten behaupten kann. 


Ein analoges, allerdings nicht zum völligen Aussterben des zahlenmäßig unterlege- 
nen Konkurrenten führendes Zusammendrängen von Allelgrenzen an Genflußbarrie- 
ren hat übrigens bereits ENDLER (1977) bei seinen auf dem Selektions-Gradienten- 
Modell beruhenden Studien beschrieben. 


Vieles spricht dafür, daß diese starke Abhängigkeit der für eine langfristige Ko- 
existenz notwendigen Mindestarealgröße vom Genflußfaktor auch dafür verantwort- 
lich ist, daß sich Organismengruppen mit extrem weit verbreiteten Diasporen (z. B. 
Moose) markant weniger in vikariierende Sippen mit historisch deutbaren Arealgren- 
zen aufgespalten haben als viele Gefäßpflanzengruppen. Besonders auffällig wird 
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dieser Unterschied beispielsweise in mediterranen Gebirgen, wo nach den Erfahrun- 
gen des Autors in Sizilien und Sardinien die Moosflora durchwegs noch wesentlich 
stärker mitteleuropäische Züge aufweist als die an Lokalentwicklungen viel reichere 
Gefäßpflanzenflora. 


Stabilität der Übergangszone gegen Fernverbreitung 


Ein prinzipieller Unterschied zwischen Selektions-Gradienten-Modellen und ,,tension 
zone“-Modellen ist, daß nur bei ersteren clinale Umweltgradienten für unterschied- 
liche Fitnesswerte der Konkurrenten und damit für die Aufrechterhaltung der 
Arealtrennung verantwortlich sind. Verpflanzungsversuche ins Areal der jeweiligen 
Konkurrenzart wären diesem Ansatz zufolge längerfristig jedenfalls zum Scheitern 
verurteilt. Komplizierter präsentieren sich die Zusammenhänge bei den „tension- 
zone“ Modellen. 
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Abb. 5: Verdrängung der von vornherein arealmäßig unterlegenen Art B durch die Art A bei gerin- 
ger Selektionsintensität (d = 0,005) und hohem Genfluß (g, = 0,1, g,, 2,, Z, g, Jeweils 0,025). Man 
beachte die autokatalytische Geschwindigkeitssteigerung dieses Vorganges in seiner Endphase. 


Da hier einseitiger Genfluß, etwa auf Grund unterschiedlicher Arealgrößen, die Art- 
frequenzen von Nachbarzellen ändern und damit dominoartige Effekte auslösen kann 
(höhere Freouenz von A > höhere Fitness von A > höherer Genfluß in nächste 
Nachbarzelle > auch hier gewinnt A die Mehrheit etc.) sind für die langfristige Co- 
existenz zweier vikariierender Arten bestimmte Mindestareale notwendig. Deren 
Größe ist dem Verhältnis Genfluß/Selektion (7) proportional. Unter dieser charak- 
teristischen Arealgröße liegende Enklaven werden auf Grund der homogenisierenden 
Kraft des Genflusses letztlich von der Nachbarart überrollt. Daß eine solcherart aus- 
reichend große Populationsinsel auf einen Schlag durch massive Fernverbreitung in- 
nerhalb des Areals eines Konkurrenten entstehen könnte, ist durch natürliche Prozes- 
se schwer vorstellbar. Am ehesten konstruierbar wäre dies allenfalls für Arten mit 
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normalerweise überaus kleinem Genflußfaktor (z. B. myrmekochore Arten), die in- 
nerhalb ihrer Areale regelmäßig größeren Populationsschwankungen ausgesetzt sind 
(z. B. durch Überschwemmungen in Augebieten). In diesem Fall wäre die zufällige 
Etablierung eines Weitwanderers an frisch zur Kolonisation freiwerdenden, von der 
Konkurrenz noch nicht erreichten Stellen denkbar. Die solcherart herangewachsene 
Inselpopulation könnte sich in der Folge - sehr niedriges Verhältnis Gen- 
fluß/Selektion vorausgesetzt — auch gegenüber einer Rückeroberung des Terrains 
durch den nun ebenfalls wieder langsam vorrückenden alteingesessenen Konkurren- 
ten als beachtlich widerstandsfähig erweisen. Möglicherweise ist das kuriose, nach 
SPETA (1982) am ehesten durch Verschwemmung erklärliche Inselvorkommen von 
Scilla bifolia L. mitten im Areal von Scilla vindobonensis SPETA in der Stopfenreu- 
ther Au so zu erklären. Hybriden zwischen den beiden chromosomal überraschend 
verschiedenen Sippen wurden bisher nicht beobachtet (GREILHUBER & SPETA 1977). 
Von besonderem Interesse wären in diesem Fall aber darüberhinaus Untersuchungen 
über den Grad der Selbststerilität sowie der Spezifität von Bestäubern. Unter der 
Voraussetzung einer (teilweisen) Selbstinkompatibilität sowie eines fehlenden Dis- 
kriminationsvermögens von Seiten der Bestäuber wäre jedenfalls die seltenere Art im 
Mischungsbereich benachteiligt, da sie mit (proportional zur Seltenheit anwachsen- 
der) Wahrscheinlichkeit durch „falschen“ Pollen bestäubt würde. Diese Situation 
hätte daher — trotz fehlender Kreuzbarkeit — überaus verblüffende Parallelen zur 
Arealgrenzbildung auf Grund einer Selektion gegen Hybriden. In beiden Fällen ist 
die häufigere Art nämlich eine „Falle“ für den Pollen der selteneren und damit ur- 
sächlich an deren (überproportionalen) Fitnessminderung beteiligt. 


Wie Abb. 5 eindrucksvoll vor Augen führt, ist in der Praxis allerdings damit zu rech- 
nen, daß die (Rück-) Eroberung eines Teilareals bei positiv frequenzkorrelierten Fit- 
nesskurven (also auch bei Hybridunverträglichkeit!) anfangs derart langsam vor sich 
gehen kann, daß dies mit unseren heutigen Methoden nicht mehr erfaßbar ist. Je mehr 
die vordringende Sippe allerdings dann die Oberhand gewinnt, umso rascher über- 
stürzen sich in der Folge die Ereignisse. Für einen menschlichen Beobachter ohne 
Kenntnis der Zusammenhänge müßten solch vergleichsweise rasche Populationszu- 
sammenbrüche nach langen Phasen unmerklicher Frequenzverschiebungen jedenfalls 
als unerklärliche Katastrophen erscheinen. 


Es ist reizvoll zu überlegen, welche Hypothesen dann für das Aussterben der letztlich 
unterlegenen Sippe aufgestellt werden könnten: Erbgutveränderungen einer der 
beiden Arten (etwa durch Mutation entstandene größere Konkurrenzkraft bei der 
vordringenden Art)? Inzuchtbedingte Degeneration bei der verschwindenden 
Minderheitsart? Eine allgemeine Klimaveränderung als Ursache für ein verändertes 
Konkurrenzgleichgewicht beider Konkurrenten? 


Übrigens kann Fitnessminderung als Folge der Unterschreitung kritischer Mindest- 
arealgrößen nicht nur bei zwei vikariierenden und dabei total konkurrierenden Sippen 
zum Aussterben eines der Kontrahenten führen. Vieles spricht dafür, daß unter sol- 


133 


chen Bedingungen auch die Wettbewerbsnachteile im steten Ringen mit einer Reihe 
anderer, jeweils nur teilweise nischenüberlappender Sippen derart gravierend werden 
können, daß sie ebenfalls mit dem Untergang lokal zu stark ausgedünnter Populatio- 
nen enden. Dieses heute leider immer aktueller werdende Naturschutzproblem wurde 
in jüngster Zeit u. a. von LAMONT & KLINKHAMER (1993) an Restpopulationen eini- 
ger vom Aussterben bedrohten Banksia-Arten Südaustraliens dokumentiert. 


Arealveränderungen bei konstanten Fitnessunterschieden sowie bei völligem 
Fehlen von Konkurrenz 


In den bisher behandelten Fällen brachten umwelt- oder frequenzbedingte Konkur- 
renzphänomene die Expansionen der beteiligten Arten letztlich dauerhaft zum Still- 
stand. Zum Vergleich sollen nun auch Modelle zur Beschreibung dauerhafter Expan- 
sionsbewegungen vorgestellt werden. In Abhängigkeit von der Konkurrenzsituation 
in den neu zu besiedelnden Räumen sind dabeı drei Fälle denkbar, welche sich in ih- 
ren Ausbreitungsgeschwindigkeiten sehr stark voneinander unterscheiden: 
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Abb. 6: Ineinanderwandern zweier überall gleich fitter Arten A und B (d = 0). Der verwendete 
mathematische Ansatz entspricht dem allgemeinen Modell (6) mit g = 0,1. 


1. Art A trifft — etwa im Zuge von nacheiszeitlichen Rückwanderungsbewegungen — 
auf Art B. Beide Sippen stehen in totaler Konkurrenz und sind überdies unter allen 
Bedingungen gleich fit (Abb. 6). Bei diesem, wohl hypothetischen Beispiel, tritt 
erwartungsgemäß langsame Vermischung ein. Je weiter diese fortschreitet, umso 
geringer wird die Ausbreitungsgeschwindigkeit beider Sippen. Der Zustand völ- 
liger Gleichverteilung wird asymptotisch erreicht. Die Verhältnisse entsprechen 
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der Situation bei der Ineinanderdiffusion zweier Gase, wobei der Genflußfaktor 
x’*g (7) einem Diffusionskoeffizienten entspricht. Am langsamsten erfolgt dieses 
Ineinanderwandern bei langlebigen Arten (K-Strategen; hier dauert der Ersatz der 
bereits etablierten Pflanzen durch Zuwanderer naturgemäß besonders lange) mit 
sehr wenig effektiver Fernverbreitung der Diasporen. Dazu gehören auch viele 
Gebirgspilanzen, u. a. die eingangs vorgestellten Sippen der Violettschwingel- 
gruppe. In solchen Fällen wird es daher im einzelnen schwer sein zu unterscheiden, 
ob die uns als stationär erscheinenden Überlappungszonen tatsächlich stabil sind 
(etwa durch positiv frequenzkorrelierte Selektion) oder ob sie nur Diffusionszonen 
mit sehr langsamen Wanderungsgeschwindigkeiten darstellen. Sind beide Sippen 
noch vollfertil miteinander kreuzbar (ohne Fitnessverlust bei den Heterozygoten) 
so läuft der Vermischungsvorgang natürlich innerhalb des gemeinsamen Genoms 
ab. Dieser Fall scheint beim Artenpaar Pulsatilla vulgaris MILL. — P. grandis 
WENDER. eingetreten zu sein, welches heute durch eine sehr weite Übergangszone 
zwischen Ostbayern und Oberösterreich verbunden ist. 


. Die aufeinandertreffenden und in totaler Konkurrenz stehenden Sippen unterschei- 
den sich — unabhängig von Populationsanteil und geographischer Position — um ei- 
nen konstanten Betrag in ihrer Fitness (Abb. 7) 
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Abb. 7: Expansion der überall fitteren Art A in das Areal der unterlegenen Art B. Dazu wurden im 
Grundmodell (6) die positiv frequenzkorrelierten Fitnessfunktionen von A und B durch Konstante 
ersetzt: W, = 1+d, W, = 1-d , d = 0,005; Simulationsannahme für den Samentransport wie bei Abb. 8 
(g, = 0,1, g, = 0,01, g, = 0,001) 


Auch solche Situationen sind nach dem sekundärem Zusammentreffen zweier vor- 
her geographisch isolierter Sippen denkbar. Beispiele hierfür wären etwa nach dem 
Abschmelzen glazialer Eisschilde ablaufende gegenläufige Rückwanderungsbe- 
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wegungen aus verschiedenen Refugien oder auch Wanderungsbewegungen der al- 
pischen und nordischen Vegetation während der Eiszeiten. Bemerkenswert ist je- 
denfalls die im Vergleich zur Expansion in konkurrenzfreie Räume (Abb. 8) über- 
raschend niedrige Wanderungsgeschwindigkeit der überlegenen Sippe. Schon 
wenn es sich um annuelle Arten handelt, welche definitionsgemäß bei jeder neuen 
Generation sämtliche Individuen ersetzen, verschiebt sich bei den angenommenen 
Parametern die Übergangszone nur um etwa 8 Zellen pro 100 Jahre. Bei lang- 
lebigen K-Strategen (wie etwa den allermeisten Gebirgspflanzen, bei denen eine 
durchschnittliche Lebenserwartung von mehr als zehn Jahren wohl die Regel ist), 
verzögert sich die Ausbreitungsgeschwindigkeit der fitteren Art natürlich noch zu- 
sätzlich. 


. Eine Art breitet sich in einen völlig konkurrenzfreien Raum hinein aus. Dieses 
Modell ist zweifellos eher eine didaktisch interessante Extremannahme als eine di- 
rekt auf die Praxis umsetzbare Simulation. Selbst der mit dem Ende der Kaltzeiten 
einsetzende Gletscherrückzug war wohl viel zu langsam, um große vegetationslee- 
re Räume zu hinterlassen. Am ehesten haben sich einige überaus konkurrenzstarke 
Neophyten an das hier simulierte Muster explosionsartiger Arealausweitung gehal- 
ten. Beispielsweise konnte man sich bei /mpatiens glandulifera ROYLE des Ein- 
druckes nicht erwehren, daß die potentiellen heimischen Konkurrenten diesem 
Neuankömmling aus dem Himalaja derart unterlegen waren, daß er an vielen eben 
erreichten Standorten tatsächlich anfangs mit exponentiellen Wachstumsraten sei- 
ne Populationen aufbauen konnte. Gerade in derart gelagerten Fällen läßt sich aber 
auch die Bedeutung der sprunghaften Fernverbreitung für derartige Ausbreitungs- 
vorgänge einer weit überlegenen Art studieren. Große Impatiens glandulifera- 
Bestände entwickelten sich nach eigener Beobachtung binnen weniger Jahre auf 
Waldschlägen fernab bereits etablierter Populationen entlang der größeren Flüsse, 
obwohl die durchschnittliche Flugweite der Schleuderfrüchte nur wenige Meter 
beträgt. Grundvoraussetzung für eine effektive Fernverbreitung sınd allerdings ne- 
ben allgemeinen Fitnessvorteilen eine gewisse Selbstfertilität (/. glandulifera ist 
nach SCHWABE & KRATOCHWIL 1991 selbstkompatibel) oder ein effektives vege- 
tatives Ausbreitungsvermögen (vgl. Veronica filiformis SMITH). Die besonders 
große Bedeutung der Fernverbreitung bei völliger Abwesenheit von Konkurrenz 
läßt sich am einfachsten dadurch erklären, daß hier jede einzelne Diaspore fernab 
vom geschlossenem Areal zu einer dauerhaften Ansiedlung (mit nachfolgend lo- 
gistischer Populationszunahme!) führt, und zwar auch dann, wenn es sich um eine 
partiell selbststerile Art mit einer sehr stark positiv frequenzkorrelierten Fitness- 
kurve handelt. 
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Abb. 8: panon einer Art in den konkurrenzfreien Raum. Nicht berücksichtigt wurde im Modell 
allerdings die Möglichkeit der Fernverbreitung, welche insbesondere bei Ruderalarten mit teilweise 
fast uneingeschränkter Selbstfertilität bzw. sehr effektiver vegetativen Vermehrung die Ausbreitung 
um Größenordnungen beschleunigt und in der Praxis häufig zur inselförmigen Gründung von neuen 
Ausbreitungszentren fernab der Stammpopulation führt. 


Modellannahme: 

Innerhalb jeder von Diasporen erreichten Zelle erfolgt logistisches Wachstum: 
[A,.,1-14,) = r*{[4,(K-[4,D/K (9) 
[A,.,) = 7*[4,]*(K-[4,VK+4,] (10) 
r = Populationswachstumsfaktor, K = Kapazität (maximale Populationsgröße) 

für die abgebildete Simulation galt: r = 2, X = 100 

Zum Zeitpunkt n heißt: 


[4,],= Frequenz von A in der ersten Nachbarzelle links, [4,], = Frequenz von A in 
der zweiten Nachbarzelle links, [A,],, = Frequenz von A in der dritten Nachbarzelle 
links. 


analog sind [4], , [A,],, , [A,],, die Frequenzen von A in den entsprechenden rechten 
Nachoarzellen. 


g, &» 8, = Genflußfaktoren für 1., 2. und 3. Nachbarzellen 
für die Simulation galt: g, = 0,1, g, = 0,01, g, = 0,001 (wie bei Abb 7) 
Die der Simulation zugrundeliegende Gleichung hat daher die Form: 


[4.1 = (0,778*((2*[4,]*(100-F4,))/100)+[4,))}+(0,1*((2*[4,1,*(100-[4,1,)/100)+[4,],)) 
+(0,01*((2*[4,],.*(100-[4,],,)/100)+ [4,],,))+(0,001 *((2*[4,],,*(100-[4,],,)/100)+ [4,],,)) 
+(0,1*((2*[4,],,*(100-[4,],,)/100)+ [4,],,))+(0,01 *((2*[A,],,*(100-[4,],,/100)+ [4,],,)) 
+(0,001*((2*[A,],,*(1 00-(4,],,)/100)+[4,],,))) 
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Selbst bei dem in unserer Simulation extrem niedrig angenommenen Populations- 
wachstumsfaktor (r = 2) erfolgt die Verschiebung der Arealgrenze (wieder unter der 
Voraussetzung der Einjährigkeit) um 98 Zellen pro 100 Jahre, also immerhin etwa 14 
mal so rasch wie in der vorangegangenen Simulation. In der Praxis wird dieser Unter- 
schied aber durchwegs noch weit größer sein, etwa wenn es sich um sehr langlebige 
K-Strategen oder (und) Arten mit hohem r handelt. Zur Veränderung des Zahlenver- 
hältnisses innerhalb einer bereits existierenden Population von leicht unterlegenen 
Konkurrenten (etwa auf Grund einer etwas geringeren Samenproduktion bei sonst 
identischen Konkurrenzeigenschaften) durch Samentransport und schließlich Selek- 
tion müssen ja zunächst einmal die etablierten Pflanzen sukzessive (und entspre- 
chende dem Verhältnis der Fitnesswerte) durch Jungpflanzen ersetzt werden. Außer- 
dem macht sıch ım Falle der Fernverbreitung einzelner Diasporen eine auf (partieller) 
Selbstinkompatibilität beruhende Fitnesssenkung bei Anwesenheit von Konkurrenten 
weit stärker bemerkbar als ım absolut konkurrenzfreien Raum. 


Die Bedeutung von Expansionsphasen für die Entstehung von Polyploiden 


Was die Evolution unserer Alpenflora betrifft, so besteht heute weitgehende Über- 
einstimmung über die entscheidende Bedeutung großräumiger Wanderungs- 
bewegungen für die Entstehung von Polyploiden aus ihren diploiden Elternsippen 
(vgl. u. a. MANTON 1937, EHRENDORFER 1949, 1962, POLATSCHEK 1966, KÜPFER 
1974, FAVARGER 1975 u.v.a.m.). Nach der heute vorherrschenden Ansicht hat dabei 
(wiederholte) geographische Isolation während der Kaltzeiten und nachfolgender 
Kontakt der sich etwas divergent entwickelten Sippen während der warmzeitlichen 
Ausbreitungsphasen zur Bildung von Allopolyploiden bzw. ,,intervarietal auto- 
polyploids“ sensu STEBBINS (1947, 1950) geführt. Da insbesonders Allopolyploide 
auf Grund der höheren Rekombinationsmöglichkeiten ihres Genoms allgemein als 
besonders adaptationsfähig gelten, sollten sie imstande gewesen sein, ihre diploiden 
Ausgangssippen teilweise zu verdrängen oder zumindest auf Grund ihrer höheren 
Anpassungsfähigkeit als erste in neu besiedelbare Räume vorzustoßen. 


Interessant wären allerdings in diesem Zusammenhang experimentelle Untersuchun- 
gen über die tatsächlichen Konkurrenzverhältnisse zwischen denjenigen diploiden 
Stammsippen, deren Verbreitung mehr oder weniger eng an eiszeitlich kaum verglet- 
scherte Refugialräume gebunden ist, und ihren polyploiden Abkömmlingen. Wären 
Polyploide nämlich prinzipiell konkurrenzkräftiger als ihre aus heutiger Sicht oft als 
„reliktisch“ eingestuften diploiden Eltern, so bleibt unklar, warum sie diese entspre- 
chend der Simulation von Abb. 7 nicht völlig aus den Refugialräumen verdrängt ha- 
ben. 
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Außerdem spricht einiges dafür, daß eine Reihe dieser Wanderungspolyploiden doch 
eher autopolyploiden Ursprungs ist. Aus unseren eigenen Untersuchungen verweisen 
wir auf Festuca norica (HACKEL) K. RiCHTER. Die hier neben der 2x-Stammform be- 
kannt gewordenen 4x- und 6x-Cytotypen unterscheiden sich allenfalls tendenziell auf 
Grund weniger quantitativer Parameter (Spaltôffnungslängen, Blattdurchmesser) von 
den Dipioiden und zeigen nicht den geringsten Hinweis auf die Beteiligung fremder 
Genome an ihrer Entstehung (PıLs 1980, DAwE 1989). Eine besonders eindrucksvolle 
Übersicht über zahlreiche ähnlich gelagerte Fälle findet sich darüberhinaus in der 
umfangreichen Arbeit von KÜPFER (1974), der übrigens explizit auf die Häufigkeit 
der Bildung von Multivalenten in der Meiose dieser polyploiden Chromosomen- 
rassen hinweist: „In many polyploids, multivalent configurations seem to be frequent, 
a fact which argues in favour of autopolyploidy and, probably, intervarietal autopo- 
lyploidy according to STEBBINS (1947, 1950).“ Letztlich muß nach SPETA (1982) 
auch die von GREILHUBER (in GREILHUBER & SPETA 1977) sowie GREILHUBER (1978) 
karyologisch recht eingehend untersuchte Scilla drunensis .SPETA als Neo- 
Autopolyploide der weiterverbreiteten S. bifolia L. gelten. Charakteristischerweise 
wagt sich erstere daher auch als einzige der österreichischen Scilla-Arten zumindest 
in E-Österreich in die größeren der inneren Alpentäler (SPETA 1987), womit sie sehr 
an die nach den Untersuchungen von PILS (1982) sowie DAWE (1989) ) wohl eben- 
falls autotetraploide „Kärnten-Steiermark-Rasse“ von Festuca pallens HOST erinnert. 


Aus all den zuletzt genannten Gründen scheint uns die Annahme einer prinzipiellen 
Überlegenheit von (Neo-) Polyploiden gegenüber ihren diploiden Stammformen beim 
derzeitigen Wissensstand nicht gerechtfertigt. Damit drängt sich unausweichlich die 
Frage auf, wie die Etablierung des ersten polyploiden Individuums in einer diploiden 
Population überhaupt vorstellbar ist. Partner für Fremdbestäubung auf der gleichen 
Ploidiestufe sind in diesem Stadium ja wohl noch nicht vorhanden, wenn man von 
der wenig wahrscheinlichen Annahme einer simultanen Entstehung gleich mehrerer 
Polyploider binnen kurzer Zeit einmal absieht. Selbstbefruchtung sollte aber die Fit- 
ness (Samenansatz) eines derartigen Neopolyploiden im ersten Stadium senken und 
seine dauerhafte Etablierung dadurch verhindern. Ein Ausweichen auf vermehrte ve- 
getative Ausbreitung zur Kompensation des mangelhaften Samenansatzes, wie dies 
etwa URBANSKA (1977a,b) bzw. URBANSKA & LANDOLT (1978) am spontan entstan- 
denen triploiden Hybriden zwischen Cardamine rivularis x C. amara feststellen 
konnten, scheint dabei eher eine interessante Kuriosität als ein allgemein gangbarer 
Weg zu sein. 
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Abb. 9: Etablierung eines Neopolyploiden im Zuge postglazialer Arealexpansion.— Zu Beginn der 
Simulation bewegt sich Art A mit den in Abb. 8 verwendeten Parametern konkurrenzfrei expandie- 
rend nach rechts. Bei Generation 41 gelingt es einer autopolyploiden Diaspore sich durch Fernver- 
breitung 18 Zellen vor der Ausbreitungsfront zu etablieren. Zum Zeitpunkt des Eintreffens der 
Hauptwelle von A hat B durch logistisches Wachstum bereits die Kapazitätsgrenze in mehreren 
Zellen rings um seine Erstansiedlungsstelle erreicht und blockiert damit unter der Annahme positiv 
frequenzkorrelierter Fitnessfunktionen (d = 0,05, g = 0,1; vor allem Hybridinkompatibilitat sollte 
dabei eine wichtige Rolle spielen!) das weitere Vordringen von A erfolgreich ab. Die weitere Ex- 
pansion in eisfrei werdendes Gebiet erfolgt von jetzt an durch den Neopolyploiden B. 


Angesichts dieser Lage erschien es reizvoll, unser Simulationsmodell einer positiv 
frequenzkorrelierten Fitness auch auf diese Situation anzuwenden (Abb. 9). Dabei ist 
nach all dem bisher Gesagtem klar, daß sich jegliche neu entstandene Mutante — so- 
ferne man ihr nicht eine überlegene Konkurrenzfähigkeit auch bei erzwungener Au- 
togamie bzw. vegetativer Vermehrung zuschreiben will — nur in extrem kleinen und 
dabei aber isolierten Populationen durchsetzen kann. Gerade dafür boten nun aber 
offenbar nacheiszeitliche Expansionsphasen die besten Voraussetzungen. Geht man 
nämlich davon aus, daß in solchen Zeiträumen auch gelegentliche Fernverbreitungen, 
etwa auf eben eisfrei gewordene Moränenrücken, Bergschultern etc. stattfinden, so 
hätten diese Populationsgründer auf Grund des Fehlens von Konkurrenten natürlıch 
Gelegenheit zu logistischem Wachstum, und zwar selbst unter der Annahme von 
„Anlaufschwierigkeiten“ auf Grund anfangs notwendiger Selbstbestaubung. 


Zum Zeitpunkt des Heranrückens der ersten Populationen der diploiden Stammart 
könnte dann dieser neopolyploide Vorposten bereits derart individuenreich geworden 
sein, daß er unter der Annahme positiv frequenzkorrelierter Fitness (in diesem Fall 
etwa auch durch Selektion gegen Hybriden!) sein inzwischen bereits etwas ausgewei- 
tetes Areal schon erfolgreich behaupten könnte. Darüberhinaus würde die weitere 
Expansion in eisfrei werdende Gebiete vom nun an der Front stehenden Neopolyploi- 
den aus erfolgen, was dessen Population natürlich exponentiell weiter ansteigen ließe 
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und damit auch gegen Unterwanderung bei hohen Genflußraten von seiten des an- 
fangs ja insgesamt zahlenmäßig weit überlegenen Diploiden unempfindlich machen 
würde. Mit diesem Modell ist demnach die in vielen Verwandtschaftskreisen nach- 
gewiesene vermehrte Präsenz von Polyploiden in eiszeitlich devastierten Gebieten 
auch ohne die Annahme einer prinzipiellen Überlegenheit von Polyploiden zwangslos 
zu erklären. Gerade hier zeigen sich aber auch die engen Grenzen derartiger Simula- 
tionen. Obwohl die allgemeine Wahrscheinlichkeit von Fernverbreitungen mög- 
licherweise sogar höher ist als heute normalerweise angenommen, sind sie doch im 
einzelnen nicht genauer vorherzusagen als Quantensprünge im Elementarteilchenzoo. 
Ohne genaueste Kenntnis von Zeitpunkt, Ort, Genom der Diaspore etc. müssen daher 
unsere Prognosemöglichkeiten ähnlich limitiert bleiben wie beim allabendlichen 
Wetterbericht. 
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Abb. 10: Simulation des Schicksals einer an wenigen Lokalstandorten fitteren Sippe im Wettbewerb 
mit einer an Normalstandorten überlegenen Konkurrenzart.- Grundvoraussetzung ist auch hier totale 
Konkurrenz (1). Weitere Simulationsannahmen: Bis zur Zelle 125 ist B (z. B. durch ein anderes 
Großklima) konstant fitter (d = 0,005), ab Zelle 126 drehen sich die Konkurrenzverhältnisse um; g, 
= 0,1. Die Simulation startet aus der durch d und g, vorgegebenen Gleichgewichtslage, d.h. A und B 
grenzen mit einer vergleichsweise schmalen Übergangszone aneinander. Lokalstandorte mit umge- 
drehten Konkurrenzverhältnissen sind die Zellen 30, 60, 90 (hier gelten innerhalb des B-Areals die 
Konkurrenzbedingungen des A-Areals) sowie 200 und 201 (hier ist B innerhalb des A-Areals die 
fittere Art). Zu Beginn der Simulation wurden diese abweichenden Lokalstandorte jeweils zu 100 % 
mit der hier fitteren Art besetzt. Wie aus dem Simulationsverlauf ersichtlich, reicht unter den ge- 
wählten Voraussetzungen eine geringe Fitnessüberlegenheit auf Grund lokalklimatischer Faktoren 
nicht aus, um großräumige Fitnessnachteile, etwa auf Grund eines die Konkurrenzart begünstigen- 
den Großklimas, auf Dauer zu kompensieren. 
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Diskussion 


Die hier vorgestellten Modelle ermöglichen eine Synthese zwischen dem Arealbe- 
griff der stark historisch orientierten, klassischen Chorologie und den Vorstellungen 
der modernen Populationsbiologie. 


Nur historisch verständlich ist zweifelsohne dıe Häufung lokal verbreiteter Arten in 
den „massifs de refuge“. Nach MERXMÜLLER (1952-54) handelt es sich dabei um 
einst wesentlich weiter verbreitete „Reliktarten“, unserer Ansicht nach um Reste ei- 
ner präglazial allgemein wesentlich kleinräumiger differenzierten Hochgebirgsflora, 
welche in den stark vergletscherten zentralen Alpenteilen stark dezimiert wurde, was 
dort zu starken Arealgewinnen der verbliebenen Sippen bzw. neopolyploider Rassen 
führte (PILS 1994). 


Durchaus aktualistischen Deutungsmodellen zugänglich sind aber diejenigen Phäno- 
mene, welche offenbar über Jahrtausende hinweg imstande waren, Arealgrenzen in 
einmal eingerasteten Positionen zu fixieren. Anstelle der historizistischen Hypothese 
einer Biotypenverarmung schlagen wir dafür ein im wesentlichen auf interspezi- 
fischer Konkurrenz beruhendes Modell vor. Dabei spricht vieles dafür, daß durch 
positiv frequenzkorrelierte Fitnessfunktionen bzw. allgemein „tension zones“ stabili- 
sierte Arealgrenzen und solche, welche durch einen Selektionsgradienten der Umwelt 
(z. B. Großklima) aufrechterhalten werden, in der Praxis oft kaum zu trennen sind. 
Lehrreich ist hier der von zoologischer Seite bereits eingehend untersuchte Fall von 
Gelbbauch- und Rotbauchunke (Bombina variegata, B. bombina). Einerseits liegen 
bei diesem durch enge Hybridzonen verbundenen, vikariierenden Artenpaar konkrete 
Hinweise für eine Selektion gegen Hybriden vor (Literaturzusammenstellung in 
SZYMURA 1993), andererseits sind beide Sıppen aber auch als indirekte Folge etwas 
unterschiedlicher Biotoppräferenzen höhenstufenmäßig geschieden. Beispielsweise 
bevorzugt die Gelbbauchunke kleinere, oft temporäre Gewässer und ist durch ihre 
dickere Haut vermutlich allgemein etwas besser an eine stärker terrestrische Lebens- 
weise angepaßt (SZYMURA 1993). 


Analoges ist natürlich auch von konkurrierenden Pflanzenarten zu erwarten. Daß sich 
„tension zones“ in Richtung auf die allgemein weniger fitte Art hin verschieben, ist 
bekannt (vgl. Abb. 7). Sie kommen letztlich dort zu Ruhe, wo der Fitnessunterschied 
zwischen beiden Sippen verschwindet (etwa auf Grund geänderter Umweltbedingun- 
gen) oder wo hinreichend effektive Genflußbarrieren dieser Verschiebung ein Ende 
setzen. Dabei muß „fitter“ keineswegs bedeuten, daß die betreffende Art in allen, 
natürlicherweise recht heterogenen Kleinstandorten (Sonnhang/Schatthang, mehr 
oder weniger basenreicher Boden....) überlegen ist. Es reicht schon, daß sie dies an 
einer Mehrzahl der Standorte (Populationen, ,,Zellen“) ist. Unter diesen Bedingungen 
machen sich nämlich an der Arealgrenze die lokalen „Fitnesshochs“ der in der Mehr- 
zahl der Zellen überlegenen Art durch einen etwas asymmetrischen Genfluß bemerk- 
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bar, was letztendlich zu deren weiterem Vorrücken führt (Abb. 10). Genau der glei- 
che Effekt würde sıch übrıgens einstellen, wenn man in einem bestimmten Simulati- 
onsbereich die Zellen mosaikartig mit wechselnden Mehrheiten von A und B besetzt. 
Auf Grund ihres ın der Summe überlegenen Genflusses würde sich auch hier die von 
vornherein häufigere Art letztlich allgemein durchsetzen. Allerdings muß für diese 
Gleichschaltung derartiger Mosaike von Kleinstandorten das Verhältnis 2. g/d 
(= Genfluß/Selektion) ein gewisses Minimum übersteigen (vgl. Abb. 5). 


Diese auch mit eigenen Sımulationen auf der Basis von Selektions-Gradienten-Mo- 
dellen ausgeloteten Gesetzmäßigkeiten sind vor allem bei der Abschätzung der Wir- 
kung großräumiger Klimagradienten auf die Pflanzenverbreitung von Interesse. An- 
gesichts der Vielzahl lokalklimatisch überaus unterschiedlicher Kleinstandorte im 
Hochgebirge mögen uns großklimatische Gradienten wie etwa die auffallend niedri- 
gen Niederschlagsmengen im nordostalpinen Refugium (vgl. PILS 1988, 1994) als 
vergleichsweise unbedeutend erscheinen. Wenn derartige großräumige Klimaver- 
schiebungen aber an der Mehrzahl der Standorte das Konkurrenzgleichgewicht von A 
(auch nur schwach!) hin zu B verschieben, reicht dies zumindest im Simulationsmo- 
dell aus um die nun unterlegene Art auch an lokalen Gunstlagen, an denen sie durch- 
aus selbst leichte Konkurrenzvorteile genießt, ebenfalls zu verdrängen (Abb. 10). 


Überaus kompliziert wird die praktische Umsetzung des hier präsentierten 
interspezifischen Konkurrenzmodells vor allem dadurch, daß die hier zwecks 
einfacherer mathematischer Behandlung vorausgesetzte absolute Nischengleichheit 
wohl ohnehin bei keinem vikarlıerenden Artenpaar tatsächlich verwirklicht ist. Daher 
sind für jede Art etwas unterschiedliche Nischenüberlappungen mit zahlreichen 
anderen Arten zu postulieren. Aus Ökologischen oder historischen Gründen können 
dabei die Areale der beiden Vikarianten recht unterschiedlich betroffen sein. Dadurch 
wird (1) natürlich in der realen Welt völlig außer Kraft gesetzt. Zellengrößen und 
damit Genflüsse sind schon aus diesem Grund nirgends innerhalb des Areals einer 
Sippe konstant, wie wir es hier der Einfachheit halber angenommen haben. Letztlich 
beschränkt sich aber unter realen Verhältnissen die Konkurrenz natürlich keineswegs 
auf systematisch nächstverwandte Vikarianten. Tatsächlich steht jede Art mit einer 
für uns derzeit noch unüberschaubaren Menge von mehr oder weniger stark 
nischenüberlappenden Sippen in mehr oder weniger heftiger Konkurrenz. Die 
Konstruktion der hier notwendigen multipolaren Konkurrenzmodelle scheitert aber 
vorläufig noch am völligen Fehlen von praktischer Feldforschung auf diesem Gebiet. 
Angesichts dieser Komplikationen wird das Hauptanwendungsgebiet derartiger 
Simulationsmodelle auf interspezifischer Ebene wohl noch eine Zeitlang vorrangig 
die prinzipielle Auslotung arealkundlicher Gesetzmäßigkeiten sein und nicht die 
detailgetreue Nachvollziehung konkreter Arealgrenzbildungen. 
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Ausblick 


Mit den hier präsentierten Computersimulationen haben wir versucht, eine Brücke 
zwischen der oft sehr historizistisch orientierten Arealkunde und den Ergebnissen der 
modernen Populationsbiologie zu schlagen. Dieser Arbeitsansatz zwingt zu mathe- 
matisch formulierbaren und damit kristallklaren Arbeitshypothesen. Die dabei ver- 
wendeten Parameter (relative Fitnesswerte, Genflußfaktoren, Grad einer eventuellen 
Hybridinkompatibilität etc.) sind Messungen zugänglich, wodurch präzisere Aussa- 
gen über die Ursachen rezenter Arealbeschränkungen möglich werden als etwa durch 
das Postulat „mangelhafter Ausbreitungstendenzen auf Grund eiszeitlicher Biotypen- 
verarmung“. Der Einsatz des Computers als „Zeitmaschine“ macht es darüberhinaus 
auch möglich, alternative Theorien über historische Art- und Arealbildungsprozesse 
auf ihre prinzipielle Sinnhaftigkeit im Lichte aktueller Erkenntnisse der Populations- 
biologie zu überprüfen. 


Was bisher allerdings fast völlig fehlt, ist die Kenntnis selbst der grundlegendsten 
Variablen solcher Modelle. Bei welchen Verwandtschaftsgruppen liegen schon ähn- 
lich konkrete Messungen der Fruchtbarkeitsminderung durch erzwungene Autogamie 
vor, wie sie von AUQUIER (1977) an Festuca durchgeführt wurden? Von welchen 
Arten wissen wir, welche durchschnittlichen Entfernungen ihre Diasporen pro Gene- 
ration zurücklegen, oder gar, mit welchen Sıppen sie in verschiedenen Teilen ihres 
Gesamtareals wie stark konkurrieren? Angesichts dieser immensen Wissenslücken 
und des bemerkenswerten Aufschwunges, welche derartige mathematisch- 
experimentell orientierte Ansätze derzeit ım angelsächsischen Raum erleben, bleibt 
jedenfalls nur zu hoffen, daß sich die botanisch-chorologischen Forschung auch hier- 
zulande bald dieses Themas annehmen wird. 


Letztlich wollen wir es aber auch nicht verabsäumen, auf die erkenntnistheoretischen 
Grenzen selbst einer derart objektivierten und mathematisierten Arealkunde hinzu- 
weisen. Sicher fehl am Platz ist die Hoffnung, daß es bald nur mehr eine Frage der 
exakten Kenntnis der Parameterwerte und einer möglichst wirklichkeitsnahen ma- 
thematischen Formelsprache sein wird, um die vergangenen und zukünftigen 10000 
Jahre Arealentwicklung der einzelnen Sippen unserer Alpenflora auf dem Bildschirm 
herunterflimmern zu lassen. Die Situation wird sich hier nie von der in der Meteoro- 
logie oder der allgemeinen Evolutionsforschung unterscheiden. Die unsichtbare 
Mauer, gegen die alle diese Wissenschaften von komplexen Systemen letztlich sto- 
ßen, wurde erstmals von Werner Heisenberg theoretisch ausgelotet. Die Konsequen- 
zen der nach ihm benannten ,,Unscharferelation“ begegnen uns auch beim Versuch 
der Simulation arealkundlicher Vorgänge auf Schritt und Tritt. Insbesondere machen 
sich mit der Zunahme der Meßgenauigkeit auch die dabei verursachten Störungen an 
den untersuchten Systemen immer drastischer bemerkbar. Übersetzt in die Alltags- 
sprache heißt dies: Unendlich genaue Messung des „status quo“ würde unendlich 
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viele Forscher mit einer unbeschränkten Speicherkapazität der beteiligten Datenver- 
arbeitung benötigen. Bekanntlich sind wir aber eben dabei, unser Weltklima schon 
mit einer vergleichsweise winzigen „Forscherbiomasse“ irreversibel zu verändern. 
Darüberhinaus machen sıch schon beı recht groben Prognoseversuchen zunehmend 
unkalkulierbare Quanteneigenschaften bei unseren Untersuchungsobjekten bemerk- 
bar. Übertragen auf das Thema unseres Artikels: Ob auf den blutverschmierten Fän- 
gen eines Adler nach seinem Kampf mit dem erbeuteten Murmeltier die Karyopse 
eines Violettschwingels gerade noch hängenbleibt, oder beim Wegfliegen doch in 
letzter Sekunde noch abgestreift wird, kann möglicherweise in Phasen der Arealex- 
pansion über die Zukunft der betreffenden Sıppe in einem ganzen Gebirgsmassiv ent- 
scheiden. 


Zusammenfassung 


Im Gegensatz zu vielen anderen Zweigen der Biologie ist in der Chorologie der Sprung von einer 
hauptsächlich deskriptiven Betrachtungsweise hin zu einem mehr experimentell-analytischen An- 


satz bisher kaum geglückt. Hauptverantwortlich dafür sind vor allem zwei Ursachen: 


1. Grundlegende Postulate der traditionellen Arealkunde sind einer experimentellen Überprüfung 
nur schwer zugänglich. Dazu gehört insbesondere die Hypothese, daß aktuelle Areal- 
beschränkungen in vielen Fällen durch „mangelhafte Ausbreitungstendenzen“ auf Grund geneti- 
scher Verarmung während Arealschrumpfungsphasen (verursacht in erster Linie durch eiszeitli- 
che Gletschervorstöße) zu erklären seien. 


2. Die Integration moderner populationsbiologischer und -genetischer Befunde in die botanische 
Arealkunde ist bisher (besonders in Mitteleuropa!) weitgehend unterblieben. 


Vor diesem Hintergrund wird versucht, ein Erklärungsmodell für allgemein anerkannte chorologi- 
sche Phänomene (etwa den besonderen Endemitenreichtum eiszeitlich unvergletschert gebliebener 
„Refugialgebiete‘) auf der Basis eines modernen populationsbiologischen Ansatzes zu erarbeiten. 
Ausgangspunkt aller Überlegungen ist die in der Populationsgenetik bereits wiederholt erprobte 
Annahme, daß aktuelle Arealgrenzen normalerweise durch zwei antagonistisch wirkende Kräfte 


fixiert werden: 


e Expansionsfördernd wirkt permanenter Genfluß (durch Diasporentransport) über die aktuellen 


Arealgrenzen hinaus. 


e Dem entgegen wirkt die Konkurrenz anderer, vikariierender Arten mit vergleichbarer (im ver- 
wendeten Modell identischer) ökologischer Nische, welche ihrerseits analoge Ausbreitungsten- 


denzen zeigen. 


Die Kombination beider Kräfte zu einer Rekursionsgleichung ermöglicht die Simulation der räumli- 
chen Entwicklung zweier total konkurrierender Arten unter den verschiedensten Annahmen und für 
viele Tausende von Generationen. 
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Ergebnisse: 


Das hier vorgestellte Grundmodell basiert auf einer positiv Frequenz-korrelierten Fitness- 
Funktion? (PFKFF) für jede der beiden konkurrierenden Arten. Unter Freilandbedingungen ist dies 
insbesondere bei obligat allogamen Arten (etwa den allermeisten Arten der Gramineengattung 
Festuca L.) zu erwarten, deren Diasporenproduktion und -vitalität der Fremdbestäubungswahr- 
scheinlichkeit und damit der Populationsdichte proportional ist. In seinen mathematischen und auch 
biologischen Konsequenzen entspricht dieses Modell den „tension zone“-Modellen der modernen 
Populationsgenetik, die stabile räumliche Grenzen zwischen Unterarten durch Fitnessreduktion auf 
Grund von Hybridunverträglichkeit erklären. Folgende allgemeingültige Schlußfolgerungen ergeben 
sich aus den durchgeführten Simulationen: 


l. 


Die Lage aktueller Arealgrenzen wird beim präsentierten Modell nicht durch aktuelle Umwelt- 
bedingungen festgelegt, sondern ausschließlich durch einen historischen Faktor, nämlich den 
Punkt des ersten Kontaktes zwischen beiden Konkurrenten (etwa nach postglazialen Expansi- 
onsphasen). 


. Wo immer PFKFFs im Spiel sind, benötigen Populationen gewisse Mindestareale, um gegen 


eine Invasion durch konkurrierende Arten stabil zu sein. Populationszusammenbrüche durch ge- 
ringfügige Unterschreitung des Minimalareals beginnen zunächst extrem langsam, zeigen aber 
in der Folge autokatalytisches Verhalten, was in der Endphase zu vergleichsweise erd- 
rutschartigen Arealeinbußen und schließlich zum Aussterben in vergleichsweise kurzen Zeit- 
räumen führt. 


Umgekehrt folgt aus Punkt 2, daß die Etablierung teilweise reproduktiv isolierter Mutanten (z. 
B. Polyploider) innerhalb von bereits existierenden Populationen mit einer PFCFF solange un- 
möglich ist, als nicht zusätzliche Effekte ins Spiel kommen, etwa eine effektive vegetative Ver- 
breitung der Mutante oder deren allgemein stark überlegene Fitness (z. B. durch Heterosis- 
Effekte). 


4. Aus Punkt 1-3 ergibt sich ein modifiziertes Erklärungsmodell für den wiederholt festgestellten 


hohen Prozentsatz von Neopolyploiden in eiszeitlich devastierten Gebirgsteilen. Eine Schlüssel- 
position nehmen dabei Gründereffekte (,,founder effects“) ein, wie sie gerade im Zuge postgla- 
zialer Expansionsphasen immer wieder zu erwarten sind. Gelegentliche Fernverbreitung einer 
mutierten (polyploiden) Diaspore in noch konkurrenzfreies Gebiet (nach dem Eisrückzug) sollte 
dort zu einer rapiden Populationszunahme durch logistisches Wachstum führen. Ist der Zeitvor- 
sprung bis zum Eintreffen der Hauptwanderungsfront des ebenfalls expandierenden Diploiden 
groß genug, so sind in der Computersimulation solche neopolyploide Gründerpopulationen be- 
reits stabil gegen eine Invasion durch den ursprünglichen Diploiden. Die weitere Expansion in 
eisfrei werdendes Gebiet erfolgt von jetzt an durch den Polyploiden, welcher solcherart selbst 
das Ausbreitungshemmnis für den vom kolonisationsfähigen Gebiet abgeschnittenen und damit 
zum „Relikt“ gewordenen Diploiden geworden ist. Hilfshypothesen wie „mangelnde Ausbrei- 


5 D. h. die relative Fitness zweier Konkurrenten A bzw. B ist ihrem jeweiligen Anteil an der Popu- 


lationssumme (A + B) direkt proportional. 
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tungsfahigkeiten“ auf Grund genetischer Verarmung werden bei einem derartigen Szenario hin- 
fällig. 

5. Zusätzlich und unabhängig vom hier vorgestellten Modell einer PFKFF können natürlich immer 
auch klimatische Gradienten zu stabilen Arealgrenzen führen. Die Bedeutung derartiger überre- 
gionaler Klimaunterschiede selbst in lokalklimatisch sehr heterogenen Gebieten wird durch eine 
abschließende, auf dem Klassischen Selektions-Gradicnten-Modcel! beruhende Simulation unter. 
strichen. Dabei zeigt sich, daß Konkurrenzvorteile einer Sippe A in eingestreuten Mikrohabita- 
ten keine Garantie dafür sind, daß sich nicht letztendlich doch die an Normalstandorten bei dem 
jeweiligen Großklima überlegene Art B überall durchsetzt. Entscheidend für den Ausgang der 
Simulation ist (wie bei Punkt 2) das Verhältnis Genfluß/Fitnessunterschied, von dem letztlich 
die Mindestgröße eines gegen Invasion durch B stabilen Mikrohabitates von A abhängt. 
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